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Abstract 

Recently a multifractal object, Q m f, was proposed to study percolation properties in a multi- 
fractal support. The area and the number of neighbors of the blocks of Q m f show a non-trivial 
behavior. The value of the probability of occupation at the percolation threshold, p c , is a func- 
tion of p, a parameter of Q m f which is related to its anisotropy. We investigate the relation 
between p c and the average number of neighbors of the blocks as well as the anisotropy of Q m f- 
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1- INTRODUCTION 



Due to the work of many physicists, and particularly to the contributions of Dietrich 
Stauffer, percolation theory has became a powerful tool in Science describing phenomena 
in many areas as geology, biology, magnetism, or social phenomena . Despite the 
enormous success of percolation it has been a theory studied in a support (lattice) that has 
a single dimension. The only references relating percolation and multifractality concern 
to the multifractal properties of some quantities of the spanning cluster at the percolation 
threshold 

n 

Recently, a model to study percolation in a multifractal was proposed p in the lit- 
erature. In fact, the authors have created an original multifractal object, Q m f, and an 
efficient way to estimate its percolation properties. In this work we study in detail the 
method to estimate p c for this multifractal and discuss the relation between p c , some 
topologic characteristics, and the anisotropy of Q m f- 

The multifractal object we develop, Q m f, is an intuitive generalization of the square 
lattice sj]. Suppose that in the construction of the square lattice we use the following 
algorithm: take a square of size 1 and cut it symmetrically with vertical and horizontal 
lines. Repeat this process n-times; at the n th step we have a regular square lattice with 
2 n x 2 n cells. The setup algorithm of Q m f is quite similar, the main difference is that we 
do not cut the square in a symmetric way. In section 2 we explain in detail this algorithm. 

The development of Q m f has a twofold motivation. Firstly, there are systems like oil 
reservoirs that show multifractal properties (f| and are good candidates to be modeled 
by such object. Secondly, there is indeed a much more general scope: we want to study 
percolation phenomena in lattices that are not regular, but that are multifractal in the 
geometrical sense. It is important to know how site percolation transition happens in 
lattices in which the cells vary in size and also in the number of neighbors. 

In this work we analyse some geometric and topologic properties of the percolation 
cluster generated on Q m f at the percolation threshold. The paper is organized as follows: 
in §2 we present the process of construction of Q m f and the algorithm for the estimation 
of Pc, in §3 we show the numerical simulations concerning the percolation threshold p c 
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and the topologic properties of Q m f'i an d finally in §4 we present our final remarks and 
comments. 

2- THE MODEL 

In this section we show the process of building of the multifractal Q m f and the key 
concepts to estimate p c . We start with a square of linear size 1 and a partition parameter 
< p < 1. For reasons that will be clear later, p = -, where s and r are integers. We call 
the multifractal built from this parameter p as the " p — Q m /\ or "(r, s) — Q m f" ■ 

The first step, n — 1, consists of two sections of the square: a vertical and an horizontal. 
Initially the square is cut in two pieces of area = and = -^- p by a vertical 
line. This process is shown in figure ^ a), where we use as an example p = ~ = |. The 
horizontal cut in which we use the same partition p is shown in figure ^b). The first 
partition of the square generates four rectangular blocks: the largest one of area (j^) 2 , 
two of area tj^w an d the smallest one of area ( y^) 2 - The difference between the largest 
area and the smallest one increases as p —> 0, moreover, the blocks get more and more 
stretched in this limit. Therefore p measures the anisotropy of Q m f- 

In the second step, n = 2, we repeat the same process of vertical and horizontal 
sections as in step 1. Generically we get 2 2n blocks after the n t/l -step. The partition 
process produces a set of blocks with a variety of areas. We call a set of all elements with 
the same area as a fc-set. At the n th -step of the algorithm the partition of the square in 
blocks follows the binomial rule: 



k=0 



The number of elements of a A-set is ('/.'• In reference Q wo *<v that as // — x all the 
Assets determine a monofractal whose dimension is dh = lim n ^ OQ l ° 9Cn s r n — . In this 

n °° log(s+r)- 

limit, the ensemble of all fc-sets engenders the multifractal object Q m f itself. 

Figure 121 shows a picture of Q m f for p = | and n = 4. We use the following code color: 
blocks of equal area have the same tonality. In other words, all the blocks of a same 
fc-set share a common gray-tone. The general view of the object shows an anisotropic, 
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heterogeneous lattice with a non-trivial topology. The anisotropy of Q m f will change the 
percolation threshold as we investigate in the next section. Before that we study in more 
detail the problem of percolation. 

The main subject of this work is the study of percolation properties of Q m f- To 
perform such a task we develop a percolation algorithm. The percolation algorithm for 
Q m f start mapping this object into the square lattice. The square lattice should be large 
enough that each line segment of Q m f coincides with a line of the lattice, this condition 
imposes that p is a rational number (it means, r and s are integers). Therefore we 
consider that the square lattice is more finely divided than Q m f- In this way all blocks 
of the multifractal are composed by a finite number of cells of the square lattice. To 
explain the percolation algorithm we suppose that Q m f construction is at step n. We 
proceed the percolation algorithm by choosing at random one among the 2 2n blocks of 
Qmf independent of its size or number of neighbors. Once a block is chosen all the cells 
in the square lattice corresponding to this block are considered as occupied. Each time a 
block of Q m f is chosen the algorithm check if the occupied cells at the underlying square 
lattice are connected in such a way to form a spanning percolation cluster. The algorithm 
to check percolation is similar to the one used in [t| 0. M. Ilf| . 

3- NUMERICAL RESULTS 

In this section we show numerical results concerning the percolation threshold. Before 
the numerics we introduce some definitions. We call lattice the square lattice underlying 



Q m f. Following the literature [l| we call p the probability of occupation of a lattice 
site. Rl is the probability that for a site occupation p there exists a contiguous cluster of 
occupied sites which crosses completely the square lattice of size L. p c is the probability of 
occupation at the percolation threshold. There are several ways j?l to define Rl- We use 
two of them: R e L is the probability that there exits a cluster crossing either the horizontal 
or the vertical direction, and R b L is the probability that there exits a cluster crossing both 
directions. At the limit of infinite lattice size R e L and R\ converge to a common value for 
the square lattice case. Besides we call p e c the value of p c estimated from R e L , it means, 
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the average p c over lattices that percolate in one direction, the horizontal or the vertical. 
And, similarly, p\ the value of p c estimated from R b L . 

Figure El illustrates the behavior of p e c and p\ for 7 different lattice sizes. Two typical 
values of p are used: p = | (solid line) and p = | (dashed line). In El (a) we show both p c 
versus 4, the upper values corresponds to p\ and the bottom values to p e c . The lattice size, 
L = (r + s) n , corresponds to 4 < n < 10, in both cases r + s = 5. The data of the upper 
branch of the figure collapse into a single curve, diversly from the lower branch. In other 
words, the curves of (2, 3) — Q m f and (4, 1) — Q m f share roughly the same pi but diverse 
p\. The reason for this feature is the strong anisotropy of (4, 1) — Q m f compared to the 
case (2,3) — Q m f. The anisotropy of the percolation cluster does not affect p c in both 
directions, but it affects p c in one direction. In fact, the percolation in both directions 
comes from an average over both directions, in such situation, any eventual anisotropy of 
the percolation cluster vanishes because of the average. 

An useful way to define p c is to take the average value p Cave = Pc ~^ Pc ■ In figure Efb) we 
plot p Cave versus ~ corresponding to the same data figure E^a). We observe in this figure 
that both curves converge to a saturation value that is not the same. The difference be- 
tween the two cases is related to the curve of p c . The anisotropy (stretching of the blocks) 
due to p implies an anisotropy in the percolation cluster. Such anisotropy determines that 
the percolation cluster does not have a correlation length independent of the direction. 
Because the percolation cluster is anisotropic the symmetry between p e c and p\ should fail. 
Therefore, we expect that for smaller p (more anisotropy in the multifractal) a greater 
difference in p Cave will appear. Figure 0] confirm this tendency for diverse values of p. 

Figure |U shows p c versus p. To obtain p c we make an average of p Cave after the saturation 
process, it means, for n > 8. The parameter p is indicated in the figure. These values are 
also shown in Table I. The straight line in the figure is the linear fitting of the data. We 
focus our attention on two features of the figure: the general tendency of decreasing p c 
with p, and the anomalous case (3, 1) — Q m f- The main tendency of decreasing p c with p 
we have discussed in connection with anisotropy. We comment the anomalous situation of 
(3, 1) — Q m f in relation to topologic properties of Q m f which are analyzed in what follows. 

The object Q m f is build of a set of blocks i with different areas, A iy and number of 
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neighbors, Q. The number Q is a central quantity in the study of topology. A topologic 
quantity, that is important in the investigation of the percolation threshold, is Cave which 
is an average of Q over Q m f- It means Cave = t 1 where iV = 2 2n is the total number 
of blocks and the sum in the numerator is performed over the full multifractal. Typically 
Cave is a number that does not increase with n, but attains a constant value after n fh 6. 
We plot in Table / this topologic quantity for n = 8 and several values of p. 

Table I 



s,r) 


(1,1) 


(4,3) 


(3,2) 


(2,1) 


(5,2) 


(3,1) 


(4,1) 


(5,1) 


Pc 


0.5929 


0.5262 


0.5262 


0.5256 


0.5252 


0.5253 


0.5243 


0.5241 


Cave 


999 


5.436 


5.436 


5.434 


5.436 


5.426 


5.436 


5.436 



We observe in the row of Cave that all these values are roughly the same, the exception 
corresponds to the case of the (3, 1) — Q m f- The conclusion we take is that the variation of 
the average number of neighbors causes the fluctuation in percolation threshold observed 
in figure EJ The p c of the (3, 1) — Q m f is a little bit greater than the other cases because 
it has a smaller average number of neighbors. In other words, because each block of the 
(3, 1) — Q m f, in the average, has less neighbors, it percolates with more difficulty than the 
tendency among its group. 

4- FINAL REMARKS 

To summarize we analyze in this work the role of the anisotropy and the average number 
of neighbors of Q m f on its percolation threshold, p c . The multifractal object is composed 
by a set of blocks with different areas and number of neighbors. As the parameter defining 
Qmf, p, goes to zero, the multifractal becomes more and more anisotropic. This anisotropy 
reflects in the percolation cluster creating an asymmetry between p e c and p c . 

The anisotropy of Q m f is evident when we compare p e c and p b c . The observed curves of 
pi show a similar behavior, in contrast to the curves of p\. Actually, the measure of p\ 
makes an average over both directions which erases any anisotropic effect of the percolation 
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cluster. This erasing effect does not exist when we measure p e c . The anisotropy of Q m f 
decreases with p (with p —> the blocks became more stretched). 

A special case in our analysis is the (3, 1) — Q m f- This case is singular compared to 
others analyzed cases. For a same number of blocks, the (3, 1)—Q m f case has less neighbors 
than the multifractal blocks corresponding to other values of p. This phenomenon is 
intrinsic to the topology of Q m f- It implies that the (3,1) — Q m f is less connected, and, as 
a consequence, it shows more difficulty to percolate. Therefore the case of the (3,1) — Q m f 
has a percolation threshold slighter greater then its neighbors in the p sequence. In a 
future work we intend to study in more detail the effect of other topologic properties on 
the percolation properties. 
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FIGURE LEGENDS 

FIG. 1: The initial step, n = 1, in the formation of Q m f. In (a) a vertical line cut the square 
in two pieces according to p. Two horizontal lines sectioning the rectangles by the same ratio 
are depicted in (b). The underlying square lattice is depicted with thick lines. 



FIG. 2: The object Q m f for p = | and n = 4, we introduce a square inside the figure to help 
the visualization. 



FIG. 3: In (a) we show p c versus for p = | (solid line) and p = \ (dashed line). The upper 
values correspond to p b c and the lower values to p e c ; it is used 4 < n < 10. In (b) we plot p ave 
versus 4- for the same data. 
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FIG. 4: The p c versus p. The chosen values of p are indicated in the figure. The straight line 
corresponds to the linear fitting. 
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This figure "figla.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0308291vl 



This figure "figlb.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0308291vl 



This figure "fig2.jpg" is available in "jpg" format from: 
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